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Comment on "Fragmented Condensate Ground 
State of Trapped Weakly Interacting Bosons in 
Two Dimensions" 

Recently Liu et al. |Q examined the lowest state of 
a weakly-interacting Bose-Einstein condensate. In addi- 
tion to other interesting results, using the method of the 
pair correlation function they questioned the valid- 
ity of the mean-field picture of the formation of vortices 
and stated that the vortices are generated at the center 
of the cloud. This is in apparent contradiction to the 
Gross-Pitaevskii approach, which predicts that the vor- 
tices successively enter the cloud from its outer parts as 
L/N (where N is the number of atoms in the trap and 
TiL is the angular momentum of the system) increases 


We have managed to reproduce the results of Rcf. 
1 1 , however a more careful analysis presented below con- 
firms the validity of the mean-field approach. 
Figures 1 and 2 show equidcnsity lines of 
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calculated within the space of N atoms and L oscilla- 
tor quanta and with the additional truncation that only 
considers single-particle components with no radial nodes 
and angular-momentum quantum numbers m = 0, . . . , 6 
(as was also assumed in Rcf. pL) Here \^l,n) de- 
notes the many-body wavefunction. In these figures 
r A = (x, y) = (3, 0) in units of the oscillator length. The 
blue colour denotes regions of lower density, while the red 
colour denotes regions of higher density. 

More specifically, in Fig. 1 we have taken TV = 40 par- 
ticles and L = 28, 32, 36, and 40, i.e., L/N = 0.7, 0.8, 0.9, 
and 1.0. Also Fig. 2 shows P(r, r A ) for L = 62, 64, 72, and 
84, i.e., for L/N = 1.55,1.6,1.8, and 2.1. These figures 
exhibit two significant differences in comparison with the 
conclusions of Ref . jjj : (i) their qualitative features are 
not independent of the radial position of the reference 
point r A , and in fact the important broken symmetry 
observed for L/N — 0.7,0.8, and 0.9 is not visible, un- 
less r A is given a significantly large value, and (ii) the 
entry of the first vortex in the region L/N < 1, as well as 
that of the second vortex in the region 1 < L/N < 1.7 is 
clearly seen to be from the outside and gradually moving 
inwards with increasing L/N. 

Finally, Fig. 1 with L /N = 1 corresponds to a single- 
vortex-like state and is slightly asymmetric with respect 
to the center because of the finite number of atoms we 
have considered (the same is also true for the slight 
asymmetries in the two-fold and the three-fold symmetric 
states shown in Fig. 2). We have confirmed numerically 
that this asymmetry decreases with increasing N. 

In summary, we disagree with the conclusions of Ref. 
jlj in the comparison of exact and mean-field description 
of the mode of entry of vortices into the condensate with 
increasing angular momentum. This does not however 
affect the conclusions of Ref. W concerning the fragmen- 
tation of a rotating condensate. 
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FIG. 1. Equidensity lines of P(r,r A ) for TV = 40 and 
L — 28, 32, 36, and 40, calculated with the reference point 
located at va = (3, 0). 
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FIG. 2. Same as Fig. 1 for N = 40 and L = 62, 64, 72, and 


84. 
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